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ABSTRACT 

We study the effects of weakly interacting massive particles (WIMP) dark matter 
(DM) on the collapse and evolution of the first stars in the Universe. Using a stellar 
evolution code, we follow the pre-Main Sequence (MS) phase of a grid of metal-free 
stars with masses in the range 5Mq < M, < BOOM© forming in the centre of a 
IO^Mq halo at z = 20. DM particles of the parent halo are accreted in the proto- 
stellar interior by adiabatic contraction and scattering/capture processes, reaching 
central densities of 0(10^^ GeV cm""^) at radii of the order of 10 AU. Energy release 
from annihilation reactions can effectively counteract the gravitational collapse, in 
agreement with results from other groups. We find this stalling phase (known as a dark 
star) is transient and lasts from 2.1 x lO^yr (M* = 600Mq) to 1.8 x lO'^yr (M* = 9Mq). 
Later in the evolution, DM scattering/capture rate becomes high enough that energy 
deposition from annihilations significantly alters the pre-MS evolution of the star in a 
way that depends on DM (i) velocity dispersion, v, (ii) density, p, (iii) elastic scattering 
cross section with baryons, ag. For our fiducial set of parameters (u, p, cro)= (10 km s"""^, 
10^^ GeV cm^'^, 10"'^^ cm^) we find that the evolution of stars of mass M„ < 4OAf0 
"freezes" on the HR diagram before reaching the ZAMS. Stars with > AOMq 
manage to ignite nuclear reactions; however, DM "burning" prolonges their lifetimes 
by a factor 2 (5) for a 600Mq {40Mq) star. 

Key words: early universe- stellar formation-dark matter 



1 INTRODUCTION 

Current observations of primordial light element abun- 
dances, baryon acoustic oscillations, distance measurements 
by means of Type la supernovae and cosmic microwave back- 
ground, all fit together to describe a Universe undergoing an 
accelerated expansion (see Spergel et al. 2007; Komatsu et 
al. 2008 and references therein). An unknown energy field, 
often referred to as Dark Energy, constitutes approximately 
75% of the total energy density, whereas the remaining is 
made of matter. However, only 15% of the latter is made 
of known particles (baryons): its majority seems to be com- 
posed of a non- visible, unknown component, commonly re- 
ferred to as Dark Matter (DM). If thermally produced in 
the hot plasma, models best fitting observations require it 
to have decoupled at temperatures much smaller than its 
mass, thus being often referred to as cold. 

In this scenario, small scale perturbations grow faster 
and detach first from the Hubble fiow, leading to a "hierar- 
chical" growth of structures starting off very small haloes in 
the young universe, and building up bigger ones by means of 



mergers. The first stars are predicted to form at a < 20 — 30 
in haloes with masses M = 10^ — 10* M©, generally re- 
ferred to as mini-haloes (see Barkana & Loeb 2001 and 
Ciardi & Ferrara 2005 for thorough reviews of the subject). 
The gas virialized in the potential wells of these systems 
has primordial chemical composition and low temperatures, 
Tvir < 10'' K; in these conditions the additional cooling nec- 
essary for the gas to collapse and form stars is provided 
by molecular hydrogen. The results of recent semi-analytic 
studies and of sophisticated 3D numerical simulations con- 
sistently indicate that the reduced cooling efficiency, to- 
gether with the absence of magnetic fields and of relevant 
angular momentum effects, inhibits gas fragmentation and 
lead to the formation of a single massive star (Omukai & 
Nishi 1998; Bromm, Coppi & Larson 1999, 2001; Abel et al. 
2000, 2002; Nakamura & Umemura 2001; Ripamonti et al. 
2002; Gao et al. 2007; O'Shea & Norman 2007; Yoshida et 
al. 2006, 2007). The mass of these first stars, often called 
Population III (Pop III) stars, is still uncertain but likely to 
be in the range 30Mq < M« < SOOMq, depending on the 
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strength of feedback effects wfiicfi control the mass growth 
by accretion on the central high-density core (Omukai & 
Palla 2003; Tan & McKee 2004). 

Since Pop III stars are predicted to form at high z or 
are hidden in the outskirts of collapsing structures at moder- 
ate z (Schneider et al. 2006; Tornatore, Schneider & Ferrara 
2007), observational evidences of their nature and proper- 
ties are still lacking. If very massive, they are thought to ex- 
plode as powerful pair-instability supernovae or to directly 
collapse to black holes after a short lifetime of a few Myr 
( |Heger fc Woosley 2002[ ). The chemical imprint they leave 
on subsequent stellar generations is difficult to identify in 
current samples of extremely metal-poor stars in the Galac- 
tic halo (Tumlinson 2006; Salvadori, Schneider & Ferrara 
2007); their signature on the reionization history is weak 
(Gnedin & Fan 2006; Choudhury & Ferrara 2006), and so 
is any feature in the low an high energy diffuse neutrino 
background (Schneider, Guetta & Ferrara 2001; locco et al. 
2005, 2008). Future probes of the nature of Pop III stars 
will come from the James Webb Space Telescope or from 
21cm telescopes which are expected to operate within the 
next decade. 

Intricate as the "standard" scenario can be, with ordi- 
nary matter only gravitationally interacting with its dark 
counterpart, there are instances where the situation could 
have been complicated by additional interactions between 
dark and ordinary matter. 

There is a flourishing zoology of models providing can- 
didates for DM particles, and we address the reader to 
Bertone, Hooper & Silk (2005) for a thorough review of 
motivations, candidates, and their properties. The currently 
favored model, which naturally complies with the require- 
ments arising from cosmological and particle physics argu- 
ments, is the lightest stable particle in a supersymmetric 
extension of the standard model of particle physics. Often 
referred to as neutralmos, these are Majorana particles cou- 
pled to baryons by means of weak interactions, with the 
most remarkable properties to be self-annihilating, and to 
have a non-vanishing scattering cross-section with standard 
model particles. 

For what concerns this paper, the consequence of these 
properties is two-fold: in environments with high enough 
density of DM particles, self-annihilation of neutralinos 
could constitute a source of energy, emitted in the form 
of radiation, which can potentially overcome the cooling 
of the gas, inhibiting or slowing down the formation of a 
proto-star. When (and if) something resembling a celestial 
object is formed, such weakly interacting massive particles 
(WIMPs) can scatter off the baryonic material and lose en- 
ergy, thus being gravitationally captured; they accumulate 
and annihilate inside the object, thus providing it with an 
(additional) energy source. 

At present, the typical DM densities are too low to pro- 
vide any dramatic, widespread effect on stellar evolution: 
recent calculations show that DM densities necessary to in- 
duce strong effects are achievable only within the central 
two parsecs of our galaxy, with similar restrictions apply- 
ing to all galaxies in the Local Universe (Fairbairn, Scott & 
Edsjo, 2008). The first star-forming mini-haloes, however, 
are smaller and denser: this provides favorable conditions 
for DM annihilation effects to play a role. 

As noticed by Spolyar, Freese & Gondolo (2008), WIMP 



annihilation in young haloes during the formation of the first 
stars could provide an amount of energy equal to the one 
dissipated by chemical cooling of the gas. Also, as noticed 
by locco (2008) and Freese, Spolyar & Aguirre (2008), the 
energy produced by DM annihilation captured inside early 
stars could even exceed the one produced by their nuclear 
burning. 

Additional effects of DM decays and annihilations have 
also been considered in the literature, with particular at- 
tention to their contribution on the reionization and ther- 
mal histories of the intergalactic medium (Mapelli & Fer- 
rara 2005; Ripamonti, Mapelli & Ferrara 2007a; Valdes et 
al. 2007) as well as on the conditions for star formation in 
the first mini-haloes (Ripamonti, Mapelli & Ferrara 2007b). 

In this paper, we aim at studying the effects of neu- 
tralino DM annihilation on stellar evolution in the early uni- 
verse. We will treat separately the effects arising from DM 
contraction during stellar collapse and from DM captured 
by scattering with the baryons: although these are clearly 
part of the same picture, the physical mechanisms are differ- 
ent and observational and experimental constraints on the 
parameters involved have different nature and reliability. 

The paper is organized as follows: in Sec. [5] we describe 
the initial conditions of the model; in Sec. |3] we introduce 
the process of adiabatic contraction and discuss the evolu- 
tion of our fiducial lOOA'/0 star in the presence of this mech- 
anism. Sec. |4] deals with the scattering/capture process and 
its impact on the evolution of the fiducial stellar model. In 
Sec. Owe discuss the dependences on the stellar mass and 
DM parameters; in Sec. [6] we summarize the effects of scat- 
tering/capture process on stellar models. Finally, in Sec. [7] 
we discuss our conclusions. We defer to the Appendixes a 
synthetic description of the stellar evolutionary code and 
include tables with relevant physical quantities of all stellar 
models considered in this study. 

Throughout the paper we work in the framework of a 
A cold DM (ACDM) cosmological model with parameters 

= 0.24, = 0.76, fis = 0.042, h = 0.73 (Spergel et 
al. 2007; Hinshaw et al. 2008) and we assume that DM is 
entirely made of neutralinos with a mass of 100 GeV. 



2 INITIAL CONDITIONS 

In this Section we describe our initial conditions for the den- 
sity profile of the dark matter halo and for the pre-Main 
Sequence (pre-MS) evolution of the star. In the rest of the 
paper we often refer to the object under investigation as a 
star, or dark star, although, in most cases, it is actually a 
proto-star. We will make the distinction clear where neces- 
sary. 



2.1 Dark matter halo 

In the present study, we implement the characteristics of 
a "standard" early star forming mini-halo, i.e. an object 
with total mass M = 10^ Mq virializing at z = 20 (see e.g. 
Abel et al. 2002; Bromm et al. 2002; Yoshida et al. 2006; 
Gao et al. 2007; Turk 2007). We assume that 82.5% of the 
total mass is DM, while the rest is baryonic; DM follows a 
standard NFW profile (Navarro, Frenk & White 1996), with 
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virial radius i?vir = -R200 = 5 x 10^" cm and concentratior0 
c = 10. We flatten this profile for radii smaller than the 
free-streaming length of DM particles (6.6 x lO^'^ cm if their 
mass is 100 GeV). 

It is worth noting that approximately IOOMq of dark 
matter, equivalent to the mass of the fiducial stellar model 
that we will discuss in the following Sections, are contained 
within a radius of ~ lO'^* cm. This qualitatively defines the 
maximum distance at which DM particles feel the gravita- 
tional pull due to the central concentration of baryons or, in 
other words, the maximum radius where adiabatic contrac- 
tion effects are relevant, as described in Sec. O 

2.2 Proto-star 

We assume that in the first mini-haloes stars form as a result 
of the collapse of metal-free gas clouds, after the cooling- 
induced fragmentation phase is completed (see e.g. Omukai 
2000 and Schneider et al. 2002). We adopt the Padova Stel- 
lar Evolution code in the version suitable for the study of 
zero metallicity stars (Marigo et al. 2001, 2003). A synthetic 
description of the code with details on the computed evolu- 
tionary tracks is presented in Appendix I All 

In order to catch the proto-star as early as possible in 
its evolution, within the convergence limits of the code and 
its physical reliability, for each stellar model we force its 
thermodynamic conditions to the tip of the Hayashi track by 
providing a density-dependent heating source; this causes an 
expansion of the proto-star and a drop of the effective tem- 
perature. We stress that this initial phase has no particu- 
lar physical meaning: any other convenient artificial heating 
source would work to this purpose. 

We prepare the initial conditions as follows. Starting 
from the configuration of a Zero Age Main Sequence (ZAMS) 
star of the same mass and adopting a primordial chemical 
compositiorl3 we then artificially expand the star towards 
the tip of the Hayashi track. During this artificial evolution 
we perform a preliminary check on the relative strengths 
of the DM and gravitational luminosities that, during this 
phase, correspond (but not exactly) to the stellar luminosity. 
For a IOOMq star, which we consider our fiducial model, 
the DM luminosity decreases as the baryonic configuration 
gets more expanded. We continue the artificial expansion 
until the ratio of DM annihilation to the total stellar one 
is Lum/L, < 0.5. At this stage, the central temperature is 
Tc ~ 5 X 10* K and the central gas density is pc ~ 10~^ g 
cm~^. The radius of the object is « lO'^* cm and, according 
to the DM profile described in the previous Section, the 
enclosed DM mass is ~ 10''^ g, only 1Q~* of the proto-stellar 
mass. It is difficult to push the proto-star beyond this point 
because of numerical problems. 

Starting from this configuration, we follow the pre-MS 
contraction phase, including the DM annihilation energy 
source term in the structure equations, as described in Sec- 
tions 13] and 131 It is important to stress that this model is 



We note that the "adiabatically contracted" DM profiles are 
almost independent of c (at least if c J; 100 — 1000). 
^ We assume a H and He mass fractions of X = 0.755 and 
Y = 0.245, respectively, according to recent BBN models 
aiocco et al. 200711 . 



physically self-consistent: if gravitational energy release is 
the only luminosity source, the model correctly reproduces 
the usual track of a "standard" , non DM-burning proto-star 
in the HR diagram. 



3 DARK MATTER CONTRACTION 

As we have discussed in the previous Section, in a proto- 
star located at the center of the halo the DM content is 
RilO"* of the total mass; as a result, the DM contribution 
to the gravitational potential is negligible, as this is largely 
dominated by the baryonic mass. Starting from the initial 
profile described in Sec. 12. II we evolve the density profile us- 
ing the so-called adiabatic contraction (AC) approximation 
(Blumenthal et al. 1986), which is based on the assumption 
that the orbital time of the particles is much longer than 
the infall time (namely, orbits never cross each other). In 
this Section, we first introduce the formalism adopted to 
implement AC in our model and then describe the evolu- 
tion of our fiducial IOOMq proto-star in the presence of DM 
contraction. 



3.1 Formalism and approximations 

The AC approximation identifies the adiabatic invariant 
M{R)R, where M{R) is the mass contained within the ra- 
dius R, as originally shown by Blumenthal et al. (1986). 
This model, which assumes spherical symmetry and circu- 
lar orbits, was improved by Gnedin et al. (2004). These au- 
thors showed that, when compared to numerical simulations, 
the Blumethal et al. model overpredicts the increase of DM 
density in the central region and that the change of the as- 
sumed invariant from M{R)R to M{R)R (where R and R 
are the current and orbit-averaged particle positions) largely 
reduces the problem. 

Gnedin et al. (2004) also estimate that, for 1Q~^ ^ 
(J?/i?vir) :S 1, 




with A ~ 0.85 ±0.05, w ~ 0.8 ±0.2. Gustaffson et al. (2006) 
confirmed these results, but showed that the values of A and 
w change from halo to halo, and that the spread is likely 
larger than the errors quoted above. However, in the fol- 
lowing we used the central values {A — 0.85, w = 0.8) from 
Gnedin et al. (2004) , as they lie well within the distribution. 

With such assumptions, the modified adiabatic invari- 
ant equation is 

-Rf[MDM,f (iif) ± Mb,f(i?f)] = iii[A/DM,i(i?i) ± Mb,i(i?i)] (2) 

where MDM,f(fi), Mb,f(-R), MDM,i(-R), and Mb,i(7?) are the 
masses of DM and baryons enclosed within a radius R, at the 
final (subscript f) and initial (subscript i) times. Given the 
initial profiles AfoM.i and Mb,i (the NFW profile described 
in Sec. 12. ip and Mb,f (the baryonic density given by the 
stellar evolution models), Eq. ([2]) can be solved iteratively 
for the radius -Rf which encloses the DM mass MDM,i(ft). 
The numerical routine which solves the equation is 
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mostly based on the public code contra by O. Gnedirlfl, 
although several adaptations and changes were necessary. 

Finally, it is important to discuss our use of Eq. 
down to ii ~ 10~^-Rvir, which is well below the limit 
(~ 10"^7?vir) where it was tested by Gnedin et al. (2004). 
Although this is definitely an untested extrapolation, we 
think our choice is well motivated and quite conservative: 
the resulting central DM density is at least a factor of 10 
lower than in the Blumenthal et al. (1986) model. It is 
worth noting that our results are in agreement with what 
recently found by Freese et al. (2008b): they study the den- 
sity profile of an AC contracted halo, adopting two differ- 
ent algorithms (based on Blumenthal's original prescription 
and on a modified method derived by Young (1980)); they 
find the two to be consistent within a factor two at baryon 
densities of C'(10~^^)g/cm'^, yielding DM densities of order 
10^''GeV/cm'^ at radii of order 10 AU. If we use the same 
initial conditions and apply the algorithm based on Gnedin's 
method to their baryonic profile, we find a DM density con- 
sistent within a factor three with the ones they obtain. 

The specific energy deposition rate due to annihilations 
of DM particles is 

dV 

where p is the local DM density, the neutralino mass, 
(av) the thermally-averaged annihilation rate. We adopt 
(av) = 3 X 10~^^ cm^ s~^ which best fits the current value 
of the DM relic density (Bertone et al. 2005). DM particles 
in the halo are strongly non-relativistic: therefore the p- 
wave term, which contributes to the annihilation rate in the 
early Universe, is negligible in astrophysical environments; 
this may lead to different values of (av) . In general, however, 
the same value we adopt is taken as "fiducial" in DM indi- 
rect search studies (see e.g. Fornengo, Fieri & Scopel 2004). 
We also assume the neutralino mass to be = 100 GeV; 
we will discuss the effects of the variation of these param- 
eters in Sec. [S] In general, only a fraction / of the energy 
released in the annihilation is emitted in form of particles 
that can be thermalized by the gas; we take / = 2/3, as for a 
typical neutralino annihilation « 1/3 of the energy goes into 
neutrinos (IBertone et al. 2005|) , and the rest of hadronic and 
electromagnetic shower induced by the primaries lastly ther- 
malizes inside the protostellar core, as estimated by Spolyar 
et al. (2008) for even lower baryonic densitiefl 



3.2 Proto— star evolution with DM contraction 

We now turn to a detailed discussion of the DM effects on 
the evolution of a 100 Mq proto-star, which we take to be 
our fiducial model deferring to Sec. [S]the study of different 
masses. We anticipate that, qualitatively, the conclusions 
and the physical picture we draw in this Section do not 
depend strongly on the assumed stellar mass. 

The initial NFW profile described in Sec. 12.11 is adia- 
batically contracted, as described in the previous Section, in 



^ See |http : //www ■ astro ■ Isa ■ umich ■ edu/~ognedln/contra/ 

* It is worth noting that the mean free path for electrons and 
photons with energy lower than the neutraUno mass m-^ , is much 
smaller than the radius of the star at any time during our analysis. 
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Figure 1. Initial NFW DM density profile of adopted M = 
IO^A/q halo (dashed line) and adiabatically contracted DM pro- 
file at the time of our initial proto-stellar phase of for the fiducial 
100 Mq star. The vertical dashed line marks the radius of the 
star at the beginning of the computation. 



order to obtain a configuration in which the baryonic com- 
ponent would correspond to our initial proto-stellar phase, 
at the tip of the Hayashi track, presented in Sec. 12.21 The 
initial NFW profile of the halo, and the AC contracted DM 
density distribution at the time we start our simulation are 
shown in Fig. [1] the dramatic enhancement (approximately 
seven orders of magnitude within the central lO^^cm) of DM 
density is the cause of the pronounced effects of DM anni- 
hilation that will be discussed in the following. 

The corresponding DM profile at this point is let evolve 
together with our stellar object, whose evolution follows a 
typical track in the HR diagram: the proto-star is totally 
convective and contracts on a Kelvin-Helmholtz time scale 
(approximately lO^yr in this phase), descending the Hayashi 
track. For reference, the track of a 100 Mq (together with 
other stellar masses) in the HR diagram is reported in Fig 

The DM annihilation luminosity (AC luminosity) be- 
comes the dominant component of the total luminosity on 
a very short timescale (~ 10 yr). Although the mass of DM 
contained inside the star at this point is only «10'^^g, 10~* 
of its baryonic mass (and 10~^ of the total mass of the halo), 
the luminosity arising from the annihilations of the DM con- 
centrated within the central ~ 2 x 10^^ cm is large enough to 
sustain the star, causing a stalling phase. This kind of object 
has been named a dark star by Spolyar et al. (2008), who 
first found that DM annihilation energy release can coun- 
teract the gravitational collapse at some point during the 
pre-stellar phase. 

Once the proto-stellar contraction has stalled, due to 
the energy released by DM annihilations, the contraction of 
the DM profile is inhibited as well. Part of the DM in the 
cusp is burned by annihilations causing a luminosity drop, 
followed by a small contraction of the baryons and DM which 
re-establishes a new equilibrium state. In reality, this se- 
quence of stable equilibrium states along the Hayashi track 
represents a continuous process which is eventually termi- 
nated when the DM density in the cusp has decreased below 
the threshold at which the energy input can no longer sus- 
tain the self-gravity pull, i.e when the annihilation timescale 
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4 DARK MATTER CAPTURE 
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Figure 2. DM density profile, truncated at the stellar radius, R,, 
at different stages during the pre— MS evolution of the fiducial 
1OOA^0 proto-star. Curves are labelled by the ratio Lqm/^*; 
conversion to corresponding time and temperatures can be read 
off Fig. [3] 



If DM is made of weakly interacting massive particles, a non- 
vanishing elastic scattering cross section between DM and 
baryons arises. WIMPs scattering off the nuclei which con- 
stitute the star lose part of their energy and some of them 
remain bound by gravitational attraction. In this Section, we 
first introduce the formalism necessary to implement scat- 
tering/capture (SC) process in the code and then we discuss 
their effects on the evolution of our IOOMq fiducial proto- 
stellar model. 



4.1 Formalism and approximations 

The capture rate, C, of DM particles by a star through scat- 
tering has been calculated by Gould (1987), and can be cast 
as follows: 



C = 47r 



where 



dRR 



dC{R) 
dV ' 



(4) 



becomes longer than the local Kelvin-Helmholtz time. The 
duration of the stalling phase, tag, is defined as the time 
needed to the AC luminosity to scale down to 50% of the 
stellar luminosity, i.e. Lom/L, — 0.5, and it is much longer 
than the typical Kelvin-Helmholtz timescale. For our fidu- 
cial set of parameters (recall this is degenerate in the ratio 
(av) /m^), the 100 Mq dark star stalls for tag = 5.3 x lO^yr. 

In Fig. [5] we show the DM density profile inside Rt at 
different stages of the proto-stellar evolution, benchmarked 
by the fraction of stellar luminosity provided by AC lumi- 
nosity, LuM/Lf The contraction and subsequent flattening 
of the central cusp with respect to the initial conditions of 
the stalling phase, where Lum/L^ = 1, together with the 
progressive shrinking and the loss of DM annihilating shells, 
is the reason of the efficiency decrease of the AC luminos- 
ity. The central enhancement of the central cusp at the 50% 
stage with respect to the 100%, dark star phase, is not able 
to compensate the loss of the external shells of annihilating 
DM. 

Fig. [3] shows the time evolution of the effective temper- 
ature for different proto-stellar masses. During the stalling 
phase, the dark star is kept stable on the HR diagram by the 
DM energy production (which is temperature independent) 
and it remains cool Tcff ~ 5 x 10^ K, as can be seen from the 
initial plateau visible in the bottom left corner of the Figure. 
When the DM cusp is exhausted, the dark star evolves along 
its pre-MS track with characteristic times longer than for 
a standard, non DM-supported, stellar model of the same 
mass, due to the additional support of AC luminosity. Fig. 
[3] restates the results that the dark star exits the stalling 
phase in tag = 5.3 x 10^ yr, and reaches the bottom of the 
Hayashi track in rnay ~ 4.4 x 10"^ yr. By comparison a stan- 
dard, non DM-supported, star of the same mass reaches the 
same point in r^^y ~ 60 yr. Table [1] shows these character- 
istic timescales for a grid of stellar mass models. 

This might have interesting implications for the radia- 
tive feedback effects of the dark star on its host and nearby 
haloes, as we will discuss in our Conclusions. 



dC{R) 
dV 



+ 



A' = 




(5) 



erf(a)], 



(Jo is the DM-baryon elastic scattering cross section, A„ 
(M„) is the atomic number (mass) of stellar nuclei, p is 
the ambient WIMP density, v is the WIMP velocity dis- 
persion, «, is the velocity of the star with respect to the 
observer, v{R) is the escape velocity at a given radius R 
inside the star, /i = m^/Mn, /i_ = (/^ — l)/2, and the sub- 
script * refers to stellar quantities. The factor rj is usually 
assumed to be •y/3/2, corresponding to the condition v* = v. 
WIMPs captured by the scattering process thermalize with 
the gas; an upper limit to the thermalization time can be 
estimated as {m^/mp){X^/vx), i.e. the WIMP-proton mass 
ratio (number of scatterings needed) times the WIMP mean 
free path divided by its dispersion velocity in the star (we 
take Vx equal to the escape velocity at the stellar surface), 
thus obtaining: 



47r 



Tth = 



R, 



,7/2 



3V2G c^o M. 



3/2 ■ 



(6) 



Assuming that the energy of the DM particles is in 
equilibrium with the gravitational field of the star, their 
number density follows a Maxwell-Boltzmann distribution 
(jGriest fc Seckel 1987|l . 



n^iR) = exp(-i? /R^), 



(7) 



where is the highest DM density achievable inside the 
star once equilibrium between capture and annihilation is 
reached after a time 
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C{av) 



(8) 



The radius within which captured DM is concentrated, once 
it has thermahzed with the star, reads 



1/2 



(9) 



where Tc and pc are the steUar core temperature and density, 
respectively. 

The energy released due to annihilations inside the star 
can be self-consistently computed once the profile n^[R) is 
known. The expression for such quantity given in Eq. ((TJ 
assumes that particles are thermahzed and equilibrium be- 
tween annihilation and capture processes has been reached. 
To take into account the transient phase before WIMPs set- 
tle to such a state, we write the annihilation luminosity as. 



LoM(t) = g(t) 47r/ / dR R n^{R){av)m^, 



(10) 



with g{t) = tanh^(f/rdyn), where rdyn=max(r;^,rth); this is a 
formal solution when t-^ > rth, and otherwise still represents 
a good approximation to deal with a transient, reducing to 
the exact solution, Ldm = Cm^f, when t > rth. Finally, / 
is the fraction of released energy absorbed within the star, 
which we take to be 2/3 for a typical neutralino annihilation, 
as discussed in Sec. [3] 

The process we have described presents some peculiari- 
ties which are worth discussing. Altough the physical energy 
source is the annihilation process, its rate is controlled by 
scattering processes which governs the capture rate. Its de- 
pendence on the background WIMP density is only linear 
(rather than quadratic, as in the case of annihilation reac- 
tions), and it depends on ctq rather than on (av). As it can 
be appreciated from Eqs. (O and (IS}, (av) and affect 
the Tx and , but within the region of the parameter space 
which is relevant to this problem they do not affect the final 
DM luminosity. 

Eq. ([5} should be integrated for each single atomic 
species in the star. However, if one relies on the current 
experimental upper limits for DM direct detection for a 100 
GeV mass neutralino, namely (To = 10~^*cm2 (Desai et al. 
20043 and Angle et al. 200®, and crj, = 4 x 10~**cm^ 
(Ahmed et al. 2004II), the capture rate is negligible for any 
species but hydrogen. In stars of primordial composition, 
such as Pop III ones, even the dependence on the coherence 
factor does not introduce any signicant contribution of 
elements other than hydrogen. Our choice for the cross sec- 
tion values is in agreement with other works in the literature, 
i. e. (|Moskalenko fc Wai 20071 IBertone fc Fairbairn 20081 
IFairbairn et al. 2008|) . 

Throughout the following we adopt « = 10 km s~^, 
which represents the virial velocity of our reference mini- 
halo with mass 10® M© at redshift z « 20. 

By integrating Eq. with a flat stellar density profile, 
one obtains a simplified expression for the capture rate. 



CocaoM.i;L— = Co'7o^ — . (H) 
m-^ Rt 

which, within the precision of experimental data, corre- 
sponds to a numerical estimate of, 

C = 9.2xlO*Vi^^ii^, 

R* JTIIOO 

having defined: 



(12) 



Pii = 



0"38 = 



10- 



mioo — 



lOOGeV ' 



and expressing M« in solar masses and R, in cm. 
It follows that, at equilibrium. 



LuM = 1.4 X 10' 



47 erg piio-38 
s R, mioo 



(13) 



(14) 



In the above expression we have taken f — 2/3 (see dis- 
cussed in Sec. (J)). Eq. (|14|l predicts that for a given mass 
M*, I/DM will grow during stellar contraction, potentially 
reaching a level able to halt the collapse. In the code, we 
have implemented a term of DM luminosity due to annihi- 
lations of captured WIMPs (SC luminosity) using Eq. (fT4ll 
multiplied by the transient factor g{t) described in Eq. (|10p . 

We emphasize that the process of DM capture is sen- 
sitive to the background DM density outside the star but 
not to the DM already accreted through adiabatic contrac- 
tion. Thus, the two processes are mutually independent. 
Moreover, scattering/capture processes can continue for long 
times: in a 10® halo, a DM luminosity of 10'*^ erg can 
be sustained for approximately 10^^ yr. 



4.2 Proto-star evolution with DM capture 

We follow the evolution of our reference 100 Mq proto-star, 
after the stalling phase due to the AC luminosity described 
in Sec. 13.21 The fiducial DM parameters are ao — lO^'^^cm'^ 
and p — lO^^GeV/cm'', which is the value set by adiabatic 
contraction in the vicinity of the star {R 10^^ cm), as 
shown in Fig. [1] It is worth noticing that this value of p 
closely matches the one predicted by 3D simulations of first 
star formation (|Turk 2007p . We will discuss the dependence 
of our model results on these parameter in Sec. [5| 

At the time of stallation, the radius of the dark star is 
R, ~ 10^'* cm. Once the DM density cusp generated through 
AC is exhausted, the energy released by DM annihilations is 
no longer sufficient to stop the gravitational collapse; also, 
the SC luminosity developed at this point is, at the equi- 
libriunjfl, approximately 10''^ erg/s, as it can be calculated 
with Eq. (|14|l . The dark star continues its evolution along 
the Hayashi line and in the pre-MS phase. 

While the star shrinks and evolves leftward on its track 
in the HR diagram, the capture rate grows as can be seen 
in Fig. m where we show the contributions to the total stel- 
lar luminosity of the different processes as a function of the 
stellar age. As shown in the Figure, despite the fact that 
SC DM annihilations are dominating the overall luminosity, 
at times >2xl0'^ yr, nuclear reactions are active producing 



a luminosity of Ln 



2 X lO''** erg/s (lO'^L©), eventually 



^ SuperKamiokande Collaboration 
6 XENON 10 Collaboration 
CDMS Collaboration 



° We note that tj/j fa 10^ yr at this stage, thus making the actual 
SC luminosity even smaller that its equilibrium value. 
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5 PARAMETER VARIATION 




0.02 



0.04 0.06 
Age [Myr] 



0.08 



Figure 3. Evolution of tlic effective temperature for a selected 
set of stellar models when only the effects of AC DM annihilations 
is considered. The dotted lines indicate the evolutionary stages of 
the stars when the AC DM luminosity has decreased to 50% and 
10% of the total stellar luminosity. 



G 



Total 

DM Adiabatic Contraction 
Gravitational Contraction 



M.= 100Mo 

DM Scattering 
Nuclear 




0.01 



0.10 
Age [Myr] 



1.00 



Figure 4. Total luminosity of our reference IOOMq star as a func- 
tion of its age. The different curves show the contribution from AC 
DM annihilations (dotted line), gravitational contraction (short 
dashed), SC DM annihilations (dot-dashed), nuclear reactions 
(dot-dot-dot-dashed). Where relevant, the processes have been 
computed for our fiducial values of ag and p (see text). 



leading the star to exhaust the hydrogen into its core and 
continue its evolution, although with longer timescales. We 
have followed the evolution of our fiducial stellar model un- 
til complete consumption of helium in the stellar core. The 
hydrogen burning lifetime is th = 4.9 Myr, to be compared 
with the Th = 2.6 Myr predicted for a Pop III star of the 
same mass in the absence of DM effects (see Table [2}. The 
helium burning lifetime remains essentially unchanged. 

Thus, the evolution of the star is slower than what ex- 
pected for a star of the same mass in the absence of DM 
capture. 

As we will show in the next Section, lower mass stars 
are more sensitive to DM effects and for this same set of 
parameters they actually stop before getting to the ZAMS. 



In this Section we will discuss the dependence of model re- 
sults on the assumed parameters. In particular, we will ex- 
plore (i) a grid of stellar masses and (ii) different DM pa- 
rameters. We will consider these in turn. 

We find that all stars stall very early in their evolution, 
when they all stand on the Hayashi track (namely, they are 
entirely convective). The evolution of the effective temper- 
ature during the AC phase for a few selected stellar models 
is shown in Fig. [3l where the dotted lines mark the evo- 
lutionary stages where the AC luminosity has decreased to 
50% and 10% of the total stellar luminosity. As can be in- 
ferred from the Figure, larger masses burn their DM content 
more rapidly than smaller ones. In Table [T] we report the 
characteristic timescale of the AC phase and the time taken 
by each star to reach the bottom of the Hayashi track with 
("THay) and without (T^ay) the contribution of AC luminos- 
ity for the whole range of stellar masses under investigation. 
The duration of the stalling phase induced by the AC lumi- 
nosity varies with stellar mass M,, ranging from 2.1 x 10"^ yr 
for a 600A/q star to 1.8 x lO'^yr for a 9Mq star. 

We have not been able to explore a wide range of values 
for {av)/m^ for problems of numerical convergence; how- 
ever, we have performed a run for our fiducial value of the 
annihilation rate and a neutralino mass = 200 GeV. Re- 



sults are reported in the Appendix IA2I and show that for 
higher values of the neutralino mass (or smaller values of the 
annihilation rate) the delaying effect of the AC luminosity is 
reduced: the smaller amount of energy per unit time makes 
the contraction of the star faster. 

Fig. [S] shows the evolution of different stellar models 
in the HR diagram. The results presented in the upper 
panel have been obtained using our fiducial DM parame- 
ters, namely ao = 10"^**cm^ and p = 10"GeV cm"^. The 
dotted lines mark the position of the star when Lum/L, — 
1 (i.e. at the beginning of the stalling phase), 0.5 and 0.1. 
As already discussed, for all the considered stellar models 
the stalling phase takes place along the Hayashi track, and 
the timing of the other benchmarked points can be inferred 
by comparison with Fig. [3] 

Following the stalling phase, these dark stars continue 
their contraction. The SC mechanism becomes more rele- 
vant as the density increases, and the SC rate grows. For 
our fiducial set of DM parameters, stars with M« J; 30Mq 
develop a SC luminosity greater than the gravitational one 
before reaching the ZAMS, and therefore do not evolve fur- 
ther on the HR diagram. This is clearly shown in the upper 
panel of Fig. O where the dashed tracks, which indicate the 
pre-MS phase, do not join the main sequence (solid lines) 
for all stellar models with Af« < 2OM0. The lower panel of 
the Figure shows the HR diagram of the same set of stel- 
lar models but assuming different DM parameters, namely 
cro=10"^^cm^ and p = 10"GeV cm"^. As expected, the 
evolution is almost unaffected by the variation of DM pa- 
rameters during the AC phase whereas the SC effects are 
significantly reduced: all the stars with M» > 5Mq, reach 
the main sequence and even the evolution of the 5Mq model 
is halted at a later time with respect to the previous case. 

In Table [5] we summarize the characteristic timescales 
which characterize the evolution of the our stellar models 
under the effect of SC assuming different DM cross sections. 
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Table 1. Characteristic times of the star relative to the phase induced by the AC luminosity. The adiabatic time tac has been defined 
as the time needed to the AC DM annihilation luminosity, Ldm to scale down to 50% of the total stellar luminosity L,. We also report 
the stellar radii Rac- The last two columns show the time required to reach the bottom of the Hayashi track from its tip with (rnay) 
and without (ri^ ) AC DM annihilation. 



M(Mq) 


TAcClO^yr) 


Rac (cm) 


THay (yr) 


-Ly(y>-) 


5 


9.7 


4.0x10^3 


3.2xl04 


1.8x10* 


7 


15 


4.6x10^3 


2.1xl04 


6.9x10^ 


9 


18 


5.0x1013 


1.8x104 


3.5x103 


12 


18 


5.6x10^3 


1.5xl04 


1.6xl03 


15 


16 


e.ixio^^ 


1.4xl04 


1.0x103 


20 


14 


6.8x1013 


1.2x104 


5.6x102 


40 


9.3 


9.0x10^3 


8.2x10^ 


2.0x10^ 


100 


5.3 


1.2x1014 


4.4x10^ 


62 


200 


3.7 


LexlOi"* 


2.5x10^ 


26 


400 


2.5 


2.1x1014 


95 


13 


600 


2.1 


2.4x1014 


45 


3.9 



as compared with the standard case where DM effects are 
neglected: once again it is evident the bigger impact of the 
SC mechanism on smaller masses and the life-prolonging 
effect of DM on all the masses. 



Fig.[6]has to be considered as a quantitative indication: 
it shows that both low and high mass stars are affected by 
DM capture. The range of stellar masses that can reach the 
ZAMS, for the fiducial DM parameter combination ((Top = 
10"^^ GeV/cm), is 40Mq < M* < IOOOMq. 



6 STELLAR MASS CONSTRAINTS 

From the above picture, it is clear that in DM-rich environ- 
ments, such as the haloes where the first episodes of star for- 
mation are expected to take place, SC luminosity may play 
a dramatic role during the early evolution of a proto-star. 
On the basis of purely quantitative arguments, this has been 
suggested by locco (2008), and Freese et al. (2008) derived 
a constraint on the mass of stars that can form under these 
conditions. Our analysis reaches conclusions that are differ- 
ent from those of the latter study as we will comment later in 
this Section. Now we want to answer the question of which 
stars will be most affected by this mechanism, and in which 
environment. Armed with the formalism developed in Sec. 
131 one may simply impose the condition Ldm^^ ^ L^ucf^t 
namely that the luminosity due to DM annihilations inside 
the star be less than the nuclear luminosity predicted at 
ZAMS for Pop III stars in the absence of dark matter. Fig. 
[S] summarizes the results of this disequation, obtained by 
applying Eq. (|14|l to a grid of stellar models at the ZAMS. 

In the region above the curve, DM luminosity exceeds 
the nuclear one and proto-stellar evolution is inhibited be- 
fore the objects reach the ZAMS, as we have discussed in 
the previous Section. Stars in this regime will "freeze" on 
the HR diagram as long as the properties of the DM dis- 
tribution around them remain the same. Below the curve 
(shaded area) stars are instead able to reach the ZAMS and 
thereafter evolve along the main sequence. 

Note, however, that the distinction between "frozen" 
and "evolving" stars has to be taken with care. In fact, as we 
have shown in Fig. [4] for our reference 100 Mq stellar model, 
the ignition of nuclear reactions can occur at very low rates, 
producing a very low nuclear luminosity, but still allowing 
the star to evolve, eventually exhausting its nuclear fuel, 
although on much longer timescales (th ~ 4.9 Myr against 
the Th = 2.6 Myr timescale predicted for a 100 Mq Pop III 
star in the absence of DM effects, see Table (2)1 . 



This is due to a change in the index k of the relation 
L* oc , following the transition of the star to a com- 
pletely adiabatic system when M, 2OOM0. In their anal- 
ysis, Freese et al. (2008) derive only an upper limit to the 
stellar mass in the presence of DM SC because they impose 
the above condition making use of the Eddington luminosity, 
which has a linear dependence on the stellar mass, instead 
of the actual nuclear luminosity at ZAMS. Therefore their 
analysis somehow under-predicts the effects of DM on low 
and intermediate mass stars. In fact, our numerical results 
show that most significant effects are obtained for objects 
with masses M, J;2OOM0, as shown in Fig. [B] This happens 
because the SC mechanism becomes efficient only when the 
stars are in a relatively advanced phase of their pre-MS evo- 
lution. 



It is also extremely interesting to notice that the mass- 
luminosity relation in low-mass stars is not particularly sen- 
sitive to metallicity. This implies that the results we have 
presented can be also applied to metal-enriched stellar pop- 
ulations, as long as they continue to form in environments 
with characteristics similar to those we have considered (see 
Sec. [21). Our results are in agreement with those of Salati & 
Silk (1989), Fairbairn et al. (2008) and Scott et al. (2007), 
who studied the main sequence evolution of WIMP-burning 
stars and observed that ZAMS stars, if "fed" with captured 
DM annihilation energy move towards the red region of the 
HR diagram, at increasing DM densities. In particular, the 
latter analysis focused on low-mass stars (A/* < 4:Mq) and 
found they eventually reach the Hayashi line for different 
DM density values at different masses (p — 10^'^GeV/cm3 
for 1 Mq), using the same current upper limit on ctq we 
adopt, and a different value for v, of the order of the rela- 
tive velocity between WIMPs and Sun, v ~ O(10^Km/s). 
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Table 2. Times needed to nuclear luminosity Lnucl to be 95% of the total L* , rz ; -tu is the time at which the hydrogen is totally exhausted 
in the core; underscript refers to the case of complete absence of any dark matter annihilation mechanism. cr3g(CT39)=10~'^*(10~^^) 
cm^, p=10"GeV/cm3. 



M(Mo) 


r|{103yr) 


rz(103yr){(739) 


rz(103yr)(<T38) 


r°H {Myv) 


T//(Myr)((T39) 


TH {Myr){a 


7 


4.9 xlO^ 


5.8 xlO-* 


Stalling 


29 


60 


Stalling 


9 


4.0 xl02 


S.lxlO'* 


Stalling 


20 


33 


Stalling 


12 


2.0x10^ 


1.6x10"' 


Stalling 


14 


18 


Stalling 


15 


1.5x10^ 


1.1x10'' 


Stalling 


11 


13 


Stalling 


20 


1.0 xl02 


6.4x103 


Stalling 


7.7 


8.6 


Stalling 


40 


41 


2.4x10^ 


2.1x10* 


4.1 


4.4 


21 


100 


22 


1.1x10^ 


4.6x10^ 


2.6 


3.0 


4.9 


200 


17 


9.6x102 


3.6x10^ 


2.2 


2.4 


3.9 


400 


14 


1.1x10^ 


3.4x10^ 


2.0 


2.0 


3.7 


600 


13 


9.7x10^ 


4.0x10^ 


2.1 


2.0 


4.1 







O 




'■-1 



5 



5.0 



4.5 

Log (Tefp, 



4.0 



3.5 



/K) 




Figure 5. The HR diagram for a grid of stellar masses. For each 
stellar model, the dashed line represents the pre-MS phase and 
the solid line represent the MS. Dotted diagonal lines mark the 
evolutionary stages when Lac/^* = li 0.5 and 0.1. The dot- 
dashed line illustrate the locus of the "freezing" points, when 
the evolution is halted by SC DM annhilation luminosity before 
the ZAMS. In the upper panel, the results have been obtained 



using our fiducial DM parameters, namely o-q 



10 



-38, 



and 



p = lO^^GeV cm~'^. The small loops in the 5 and 7 Mq models 
are due to the effect of WIMPs thermalization, which results in 
an effective delay of the effects of SC DM annihilation luminosity. 
In the lower panel, the stellar models have been run assuming the 
same DM density but (jg = lO^^^cm^. 




10 



100 



M./Me 



Figure 6. Dark matter constraints on stellar mass. The shaded 
area represents the region of the parameter space where stars 
can reach the ZAMS and evolve. The vertical axis on the left 
represents the quantity crop and on the right the corresponding 
values of p if ctq = lO^^^cm^. In the region above the curve, 
stars are prevented from reaching ZAMS and freeze on the HR 
diagram, as it is shown in Fig. \5\ 



7 CONCLUSIONS 

We have studied the effects of WIMP dark matter annihi- 
lation on the first stars in the Universe. As initial condition 
of our model, we consider a dark matter halo with mass 
10^ Mq at 2 = 20, i.e. a typical mini-halo where the first star 
formation episodes are expected to occur. We have treated 
separately the mechanism of adiabatic contraction (AC) and 
scattering/capture (SC) and highlighted their effects on the 
pre-Main Sequence (pre-MS) phase of stellar objects with 
masses 5Mq < M^, < 6OOM0 formed from the collapse of 
metal-free gas clouds. We find that: 

• Early in the proto-stellar evolution, the luminosity pro- 
duced by DM annhilations during the AC regime induces 
a transient stalling phase: all the proto-stars become dark 
stars for characteristic times ranging from 2.1 x lO^yr for a 
6OOM0 star to 1.8 x Iff^yr for a 9Mq star. 

• The stalling phase occurs when the stars are on the 
Hayashi track. The AC luminosity does moderate the ef- 
fective temperature of the star, by enabling an equilib- 
rium state at the early evolutionary stages, characterized 
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by larger radii, when the effective temperature is «5xl0'' 
K. 

Later in the evolution, the capture of WIMPs by means 
of scattering on the baryonic matter of the star becomes high 
enough that WIMPs accumulated and annihilating inside it 
produce enough energy to keep the object at equilibrium. 
The details of this depend on environmental conditions of 
DM, elastic scattering cross section between WIMPs and 
baryons, and mass of the star. 

For our fiducial set of parameters {v=10 km/s, 
p=10"GeV/cm^ cro^lQ-^^^cm^), we find that: 

• Small and intermediate mass stars (M* < 4OM0) are 
most affected by SC luminosity and their evolution is halted 
on the HR diagram before reaching the ZAMS. 

• Dark stars can be supported by SC luminosity as long 
as the environmental conditions remain unaltered. This is 
because, unlike the AC mechanism, SC luminosity depends 
on the flux of DM particles streaming through the star from 
outside, thus drawing from a virtually unexhaustable reser- 
voir: dark stars remain "frozen" on the HR diagram. 

• Stars with masses > 4OA/0 manage to ignite nuclear re- 
actions, and go through the Main Sequence supported by an 
additional energy source: dark matter "burning" prolonges 
their lifetimes from a factor 2 for a 6OOM0 to a factor 5 for 
a 4OM0 star. 

These conclusions depend on the assumed dark matter 
parameters and on the specific environment where the first 
stars are expected to form. However, they do not strongly 
depend on the assumed primordial chemical composition of 
the stars. Thus, they can be applied to more evolved stellar 
populations as long as the characteristics of the environ- 
ment where these stars form are similar to the ones we have 
considered. 

The existence of dark stars can have many interesting 
consequences for a number of issues. In fact, once they are 
frozen on the HR diagram, they have effective temperatures 
in the range 10* — 10^ K and provide a continuous source 
of UV photons. This could have interesting consequences 
for the radiative feedback on the parent and neighbouring 
haloes as well as on the reionization history. The duration of 
this dark stellar phase is difficult to estimate, as it depends 
on the persistence of high dark matter densities around the 
stars. Therefore, their fate is strictly related to the evolution 
of their parent dark matter haloes and their merger histories. 

Of course, the present analysis represents only a first 
step of more refined future studies; however, it opens a rel- 
atively novel window on high redshift star formation. Pro- 
gresses on the issues discussed here might lead to consider- 
able understanding of the signatures that the yet mysterious 
dark matter particles have unmistakably left on the stars 
that formed in the baby Universe. 
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APPENDIX A: STELLAR EVOLUTION 
Al The code 

The modifications introduced in the Padova stellar code to 
follow the evolution of zero-metallicity stars are fully de- 
scribed in Marigo et al. (2001, 2003) to which the reader is 
referred for details. 

In the following we summarize the main input physics 
and further implementation required to follow the pre-main 
sequence phase. 

Radiative opacities are from the OPAL group, Rogers 
& Iglesias (1992) and Iglesias & Rogers (1993), for temper- 
atures T >10*K, and from Alexander & Ferguson (1994) 
for T <10'*K. Conductive opacities of electron-degenerate 
matter are from Itoh et al.(1983). 

The equation of state for temperatures higher than 
10^ K is that of a fully-ionized gas. At high densities. 
Coulomb interactions are introduced adopting the prescrip- 
tion by Straniero (1988) 

Nuclear reaction rates are from the compilation of 
Caughlan & Fowler (1988) while, energy losses by pair, 
plasma, and bremsstrahlung neutrinos, are from Haft, Raf- 
feh & Weiss (1995). 

The energy transport in the outer convection zone is 
described according to the mixing-length theory of Bohm- 
Vitense (1958) with a mixing length parameter, a = 1.68, 
obtained from the calibration of the solar model by Girardi 
et al. (2000). 



The extension of convective boundaries is estimated 
with the standard Schwarzschild criterion. Adopting more 
sophisticated schemes such as semiconvection and/or con- 
vective overshoot (e.g. Bressan et al. 1981 and Alongi et al 
1991) would not change the results. 

Effects of stellar rotation and/or magnetic fields have 
not been considered in this exploratory work. 

We finally remark that for zero-metallicity stars the 
abundance equations need to be solved simultaneously, for 
both the H and He burning reactions without any assump- 
tion for nuclear equilibria. This is performed with a semi- 
implicit extrapolation scheme (Bader & Deuflhard 1983). 

A2 Evolutionary Tracks 

We have computed evolutionary tracks for initial masses in 
the range 5 Mq to 600 Mq, starting from the pre-main se- 
quence phase and covering the central Hydrogen and Helium 
burning phases. The adopted initial composition consists of 
a mixture of just hydrogen and helium, resembling the lack 
of metals in the early Universe ppcco et al. 2007|) . with mass 
fractions of X — 0.755 and Y = 0.245, respectively. 

The evolution is performed at constant mass, i.e. ne- 
glecting both mass accretion during the pre-main sequence 
phase and mass loss by stellar winds in the later stages. 

Five different sets of tracks have been considered, de- 
pending on the assumptions concerning the DM parameters. 
The results are summarized in Tables Al to A5. For selected 
stages during the evolution we report: age, surface luminos- 
ity L, effective temperature T^ff, radius i?*, central values 
of density pc, temperature Tc and hydrogen fraction Xc, 
and the fractional luminosity provided by nuclear reactions 

(LaTuc/L*)- 

Table Al refers to stellar evolution without DM effects 
(the standard case). In Table A2 we consider only the ef- 
fects of annihilation of DM in the adiabatic contraction 
phase. Tracks in tables A3 and A4 also take into account 
the additional role of DM scattering/capture processes for 
two different values of the parameter ao, ao = W~^^cm^ and 
ao — 10~^*cm'^. Finally in Table A5 we have analysed the 
influence of a different neutralino mass (200 GeV instead of 
100 GeV) on the adiabatic contraction phase. 

The different stages selected in the Tables have the fol- 
lowing meaning. 

The Starting Point marks the beginning of the evolu- 
tionary sequences. In the standard case this is the first point 
where the model is fully and consistently supported by the 
gravitational energy release, while in the AC cases it is the 
model where the total luminosity of the star is equalized by 
the AC luminosity of its current baryonic configuration. At 
these luminosities the contraction timescale are of the order 
of a few years, so that some differences in the starting lumi- 
nosity do not affect the subsequent evolutionary timescales. 

The Hayashi Minimum Luminosity Model corresponds 
to the stage of minimum luminosity during the descent along 
the Hayashi locus. Soon after this point the proto-star moves 
towards the main sequence at nearly constant luminosity. 
This is a convenient point for comparing the contraction 
timescales of the different sets and their sensitivity to DM 
adiabatic contraction. 

The points labelled Lac /L* =50% and 10% indicate the 
stages when the AC luminosity contributes to the 50% and 
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10%, respectively, of the total stellar luminosity. The com- 
plementary energy source comes from gravitational contrac- 
tion. When Lac/L*=10%, the adiabatic contraction phase 
is essentially over. 

The stage Ljvuc./L,= 95% indicate the point where the 
nuclear energy source provides 95% of the stellar luminos- 
ity. In absence of DM scattering/capture this point indicates 
the beginning of the major Hydrogen nuclear burning phase. 
On the contrary when DM scattering/capture becomes ef- 
ficient, this point may be reached at an advanced stage of 
nuclear burning, or even not reached at all, for stars that 
suffer SC stalling. The contribution of the SC energy source 
makes the nuclear burning to occur at a slower rate, thus 
prolonging the lifetime of the star by a significant amount. In 
some circumstances the star is practically totally sustained 
by the SC luminosity with a negligible nuclear burning. In 
all these cases we provide the central hydrogen fraction and 
the fractional nuclear luminosity. These quantities may help 
the reader to evaluate the importance of the effect and, for 
the stalled stars, to roughly estimate the duration of the 
phase. 

Finally, the last stage in the tables corresponds either to 
central Hydrogen exhaustion, or to the stalled model at an 
age of 5 Myr. We remark that, in the latter case, the stellar 
tracks have been evolved for a much longer time (larger than 
the standard Hydrogen burning lifetime) to properly check 
the stalling condition. In Tables A3 and A4 we do not report 
the first stages because they are identical to those in Table 
A2. 

A few remarks on the effects of the AC and SC mecha- 
nisms are worthy at this point. 

The AC term is important in the early contraction 
pre-main sequence phase. By comparing AC and Standard 
models we notice that, while the absolute duration of this 
phase increases as the mass decreases, the relative effect (i.e. 
the ratio between lifetimes at the Hayashi Minimum for a 
given mass) increases strongly with the stellar mass. In any 
case the AC phase does not particularly affect the total nu- 
clear burning lifetime. Actually some Standard models show 
slightly longer Hydrogen burning lifetimes compared to the 
corresponding AC cases. This reflects the fact that the be- 
ginning of the main sequence is reached with small struc- 
tural differences (e.g. location of convective boundaries), 
that propagates during the subsequent evolution. 

The SC term is never important in the pre-main 
sequence contraction phase since the large stellar radii 
strongly limit the SC luminosity (directly via Eq. (|14|) and 
indirectly through the corrective term due to thermalization 
time, Eq. (O). 

This paper has been typeset from a TjrjX/ I^TgX file prepared 
by the author. 
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Table Al. The values for "standard", metal-free stars evolving without DM annihilation. 



M/Mq 


Age(Myr) 


Log(L/LQ) 


LogfTg//) 


R*(cm) 


Log(pc) 


Log(To) 








starting P. 


aint 




















5 


0.0000000 


4.762 


3.646 


2.8556 


10l3 


-6.317 


4.898 


- 


- 




7 


0.0000000 


4.854 


3.651 


3.0992 




-6.269 


5.008 


- 


- 




9 


0.0000000 


5.096 


3.649 


4.1400 


10l3 


-6.611 


4.962 


- 


- 




12 


0.0000000 


5.223 


3.652 


4.7252 


10l3 


-6.650 


5.014 


- 


- 




15 


0.0000000 


5.305 


3.656 


5.0942 


ioi3 


-6.624 


5.068 


- 


- 




20 


0.0000000 


5.416 


3.660 




10^^ 


-6.604 


5.129 


- 


- 




40 


0.0000000 


5.769 


3.669 


8.1776 


10l3 


-6.692 


5.211 


- 


- 






n nnnnnnn 
u.uuuuuuu 






1.1439 


ioi4 












200 


0.0000000 


6.471 


3.715 


1.4849 


ioi4 


-6.495 


5.466 








dOO 


0.0000000 


6,7^2 
T.inti 


3,73G 


1,9331 


loi-i 

11,11 


-6, 119 
171 


.^.,,'17 








Hayashi M 


inimum Lumi 


nosity Model 


















5 


0.0180837 


2.643 


3.739 


1.6218 


10l2 


-2.038 


6.237 








7 


0.0069470 


3.103 


3.730 


2.8721 


10l2 


-2.640 


6.129 








9 


0.0035714 


3.432 


3.724 


4.3114 


ioi2 


-3.058 


6.058 








12 


0.0016542 


3.790 


3.717 


6.7326 


10l2 


-3.558 


5.971 








15 


0.0010506 


4.057 


3.713 


9.3245 


10l2 


-3.838 


5.931 








20 


0.0005610 


4.375 


3.708 


1.3782 


10l3 


-4.261 


5.861 








40 


0.0002053 


5.034 


3.703 


3.0049 


10l3 


-4.917 


5.776 








100 


0.0000618 


5.742 


3.707 


6.6582 


10l3 


-5.633 


5.672 








200 


0.0000261 


6.265 


3.726 


1.1135 


lo" 


-5.990 


5.632 








400 


0.0000132 


6.748 


3.754 


1 . 7054 


lo" 


-6.154 


5.644 








600 


0,0000039 


6,997 


3,7.'7 


2,2,'>03 


lo" 


-6.372 


.5,608 








L., - !).->,-; 


5 


0.7510101 


2.925 


4.453 


8.3722 


lolo 


2.090 


7.676 


0.75194 


9.5000 10~ 


1 


7 


0.4728800 


3.378 


4.552 


8.9459 


lolo 


2.114 


7.777 


0.75154 


9.5001 10- 


1 


9 


0.3380348 


3.705 


4.624 


9.3718 


lolo 


2.139 


7.851 


0.75111 


9.5000 10- 


■ 1 


12 


0.2012758 


4.072 


4.703 


9.9462 


lolo 


2.169 


7.931 


0.75197 


9.5000 10" 


1 


15 


0.1537508 


4.340 


4.762 


1.0292 


loii 


2.200 


7.992 


0.75210 


9.5000 10- 


■ 1 


20 


0.1005507 


4.699 


4.829 


1.1443 


lo" 


2.201 


8.051 


0.75305 


9.5000 10- 


1 


40 


0.0407115 


5.456 


4.921 


1.7910 


lo" 


1.998 


8.100 


0.75440 


9.5001 10~ 


1 


100 


0.0220905 


6.161 


4.988 


2.9595 


lo" 


1.748 


8.138 


0.75480 


9.5002 10- 


■ 1 


200 


0.0167265 


6.599 


5.018 


4.2737 


loii 


1.585 


8.159 


0.75480 


9.5224 10" 


1 


400 


0.0136819 


6.995 


5.033 


6.2788 


lo" 


1.443 


8.177 


0.75490 


9.5000 10~ 


■ 1 


600 


0.0130140 


7.208 


5.039 


7.8188 


loii 


1.359 


8.185 


0.75490 


9.5001 10- 


1 


Central Hydrogens 


5 


55.5662916 


3.336 


4.504 


1.0655 


lo" 


3.214 


7.974 


0.00000 


0.0000 




7 


29.3981445 


3.728 


4.567 


1.2473 


lo" 


2.916 


8.039 


0.00000 


0.0000 




9 


19.8050932 


4.061 


4.602 


1.5618 


loii 


2.758 


8.077 


0.00000 


0.0000 




12 


14.0895636 


4.456 


4.638 


2.0846 


lo" 


2.604 


8.110 


0.00000 


0.0000 




15 


10.5241264 


4.716 


4.666 


2.4681 


lo" 


2.520 


8.129 


0.00000 


0.0000 




20 


7.6663778 


5.034 


4.695 


3.1161 


loii 


2.411 


8.150 


0.00000 


0.0000 




40 


4.1366357 


5.666 


4.740 


5.2562 


lo" 


2.216 


8.188 


0.00000 


0.0000 




100 


2.6233011 


6.311 


4.748 


1.0603 


10l2 


1.981 


8.219 


0.00000 


0.0000 




200 


2.1888969 


6.718 


4.736 


1.7944 


10l2 


1.824 


8.236 


0.00000 


0.0000 




400 


1.9639194 


7.088 


4.659 


3.9227 


10l2 


1.693 


8.254 


0.00000 


0.0000 




600 


2.0630663 


7.418 


4.785 


3.2025 


10l2 


1.594 


8.268 


0.00000 


0.0000 





Table A2. AC mechanism only 



M/Mq 


Age(Myr) 


Log(L/LQ) 


LogfTg//) 


R*(cm) 


Log(pc) 


Log(To) 












starting P. 


aint 
























5 


0.0000000 


5.022 


3.637 


4.0174 


10l3 


-6.924 


4.712 




- 




- 




7 


0.0000000 


5.151 


3.642 


4.5638 




-6.930 


4.798 




- 




- 




9 


0.0000000 


5.247 


3.645 


5.0209 


10l3 


-6.943 


4.860 




- 




- 




12 


0.0000000 


5.357 


3.648 


5.6097 


10l3 


-6.945 


4.922 




- 




- 




15 


0.0000000 


5.444 


3.651 


6.1288 


10l3 


-6.941 


4.968 




- 




- 




20 


0.0000000 


5.556 


3.655 




ioi3 


-6.919 


5.027 




- 




- 




40 


0.0000000 


5.837 


3.667 


8.9564 


10l3 


-6.845 


5.161 




- 




- 






n nnnnnnn 
u.uuuuuuu 






1.2491 


ioi4 
















200 


0.0000000 


6.547 


3.715 


1.6202 


ioi4 


-6.622 


5.424 












dOO 




O.Rll 


:^,733 


2,0096 


loi-i 


-r.,7 76 


















7,1121 






11)11 


















50% 
























5 


0.0096533 


3.074 


3.717 


2.9553 


10l2 


-3.216 


5.937 












7 


0.0151310 


3.201 


3.722 


3.3305 


10l2 


-3.086 


6.018 












9 


0.0177185 


3.449 


3.724 


4.3954 


ioi2 


-3.043 


6.067 












12 


0.0175911 


3.886 


3.726 


7.2124 


10l2 


-3.174 


6.095 












15 


0.0163967 


4.187 


3.728 


1.0079 


10l3 


-3.375 


6.086 












20 


0.0144465 


4.546 


3.748 


1.3908 


10l3 


-3.608 


6.081 












40 


0.0092501 


5.259 


3.766 


2.9110 


10l3 


-4.295 


5.992 












100 


0.0053198 


6.040 


3.943 


3.1780 


10l3 


-4.055 


6.204 












200 


0.0036779 


6.514 


4.041 


3.4862 


lo" 


-3.914 


6.326 












400 


0.0024905 


6.935 


4.094 


4.4213 


10l3 


-3.871 


6.406 












600 


0,0020670 


7,164 


4,110 


,5,3619 


ioi3 


-3.881 


6,439 














5 


0.0518843 


2.801 


3.763 


1.7432 


10l2 


-1.456 


6.419 












7 


0.0419409 


3.354 


3.891 




jq12 




6.559 












9 


0.0341297 


3.719 


3.969 


1 9446 


j^q12 


"l^71 














12 


0.0287672 


4.110 


4.070 


1.9156 


10l2 


-1.305 


6.726 




- 




- 




15 


0.0250591 


4.388 


4.134 


1.9676 


ioi2 


-1.293 


6.788 












20 


0.0215014 


4.719 


4.221 


1.9227 


10l2 


-1.205 


6.888 












40 


0.0136990 


5.392 


4.365 


2.1521 


10l2 


-1.184 


7.036 












100 


0.0076639 


6.099 


4.454 


3.2231 


10l2 


-1.376 


7.098 












200 


0.0057030 


6.553 


4.527 


3.8831 


10l2 


-1.310 


7.194 












400 


0.0042986 


6.959 


4.571 


5.0526 


10l2 


-1.316 


7.258 












600 


0.0037723 


7.180 


4.586 


6.0964 


10l2 


-1.352 


7.281 












Hayashi M 


inimum Lumii 


losity Model 






















5 


0.0318454 


2.650 


3.739 


1.6363 


10l2 


-2.062 


6.232 












7 


0.0207489 


3.115 


3.730 


2.9142 


10l2 


-2.650 


6.130 












9 


0.0177888 


3.449 


3.724 


4.3936 


10l2 


-3.036 


6.069 












12 


0.0153947 


3.814 


3.717 


6.9172 


ioi= 


-3.536 


5.982 












15 


0.0144330 


4.089 


3.713 


9.6710 


10l2 


-3.842 


5.936 












20 


0.0124400 


4.410 


3.707 


1.4378 


loi-^ 


-4.283 


5.859 












40 


0.0082365 


5.094 


3.699 


3.2837 


10l3 


-5.150 


5.707 












100 


0.0043604 


5.805 


3.704 


7.2797 


10l3 


-5.820 


5.613 












200 


0.0024627 


6.297 


3.721 


1.1845 


ioi4 


-6.130 


5.587 












dOO 


1), 0000.12^ 


0,71^ 


:V730 


1,^37,7 


lo" 


-r>,300 
















(J.ltlMt 1 








11,11 


















= 95% 
























5 


0.6821600 


2.928 


4.452 


8.4617 


lolo 


2.079 


7.672 





.75305 


9 


.5000 10~ 


1 


7 


0.4372555 


3.380 


4.551 


9.0212 


lolo 


2.105 


7.774 





.75265 


9 


.5000 10~ 


1 


9 


0.3242375 


3.707 


4.623 


9.4319 


lolo 


2.133 


7.849 





.75211 


9 


.5000 10~ 


1 


12 


0.2193152 


4.071 


4.703 


9.9367 


lolo 


2.169 


7.931 





.75208 


9 


.5000 10~ 


1 


15 


0.1737158 


4.342 


4.762 


1.0320 


loii 


2.198 


7.992 





.75206 


9 


.5000 10^ 


1 


20 


0.1152623 


4.699 


4.829 


1.1440 


loii 


2.201 


8.051 


0, 


.75310 


9 


.5000 10~ 


1 


40 


0.0510427 


5.455 


4.921 


1.7890 


loii 


1.999 


8.100 





.75440 


9 


.5001 10~ 


1 


100 


0.0278809 


6.161 


4.988 


2.9578 


loii 


1.749 


8.138 





.75480 


9 


.5001 10~ 


1 


200 


0.0208053 


6.599 


5.018 


4.2722 


loii 


1.586 


8.159 





75480 


9 


.5001 10~ 


1 


400 


0.0168216 


6.995 


5.033 


6.2819 


loii 


1.442 


8.177 





75490 


9 


.5001 10~ 


1 


600 


0.0152766 


7.208 


5.039 


7.8170 


loii 


1.360 


8.185 


0, 


.75490 


9 


.5080 10~ 


1 


Central Hydrogen= or 


Scattering Stallation model 




















5 


56.1465439 


3.340 


4.505 


1.0660 


lo" 


3.196 


7.983 





.00000 




0.0000 




7 


29.1142533 


3.725 


4.569 


1.2323 


loii 


2.929 


8.037 


0, 


.00000 




0.0000 




9 


19,2892971 


4.046 


4.607 


1.4988 


loii 


2.772 


8.074 





00000 




0.0000 




12 


13,2154438 


4.420 


4.646 


1.9311 


loii 


2.624 


8.107 





00000 




0.0000 




15 


10.7956074 


4.734 


4.663 


2.5626 


loii 


2.515 


8.131 


0, 


.00000 




0.0000 




20 


7.6619203 


5.035 


4.695 


3.1251 


loii 


2.415 


8.151 





.00000 




0.0000 




40 


4.1429360 


5.667 


4.740 


5.2378 


lo" 


2.219 


8.189 


0, 


.00000 




0.0000 




100 


2.6907600 


6.325 


4.740 


1.1194 


10l2 


1.986 


8.222 





00000 




0.0000 




200 


2.1910250 


6.720 


4.712 


2.0064 


10l2 


1.831 


8.239 





00000 




0.0000 




400 


2.2893942 


7.330 


5.117 


6.2717 


loii 


1.663 


8.269 





.00000 




0.0000 




600 


1.8355234 


7.289 


4.543 


8.4022 


10l2 


1.612 


8.261 





.00000 




0.0000 
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Table A3. AC and SC mt-cliauisins active: f)^= lO^i OoV/rm'^ <t()= m-'-^''cm'^ 



M/Mq 


Age(Myr) 


Log(L/LQ) 




R* (cm) 




Log(To 


) Xe 






Hayashi M 


inimum Luminosity Model 
















5 


0.0321649 


2.650 


3.739 


1.6353 10^^ 


-2.056 


6.233 


_ 


_ 




7 


0.0207489 


3.115 


3.730 


2,9142 10^2 


-2.650 


6,130 


_ 


_ 




9 


0.0177888 


3.449 


3.724 


4,3936 10^2 


-3.036 


6,069 


- 


- 




12 


0.0153947 


3.814 


3.717 


6.9172 10^^ 


-3.536 


5.982 


- 


- 




15 


0.0144330 


4.089 


3.713 


9.6710 10^^ 


-3.842 


5.936 


- 


- 




20 


0.0124400 


4.410 


3.707 


1.4378 10^^ 


-4.283 


5.859 


- 


- 




40 


0.0082365 




3.699 


3,2837 10^^ 


-5.150 


5,707 


- 


- 




100 


0.0043604 


5.805 


3.704 


7,2797 10^'^ 


-5.820 


5,613 


- 


- 






0.0024627 


6.297 




1,1845 10^^ 












300 


0.0016023 


6.565 


3.732 


1.5320 10^^ 


-6.254 


5.585 


- 


- 




400 


0.0009528 


6.748 


3.739 


1.8375 Ifll* 


-6.366 


5.574 








600 


0.0004468 


6.998 


3.756 


2.2582 10^^ 


-6.386 


5.603 










= 95% 


















7 


58.2469694 


3.692 


4.559 


1,2448 10^1 


2.532 


7,974 


0,11475 


9,5000 10" 


1 


9 


30.9662677 


3.992 


4.607 


1.4085 10^^ 


2.378 


7.985 


0.15375 


9.5000 10~ 


1 


12 


15.8442989 


4.330 


4.663 


1.6080 10^^ 


2.226 


7.997 


0.21329 


9.5000 10~ 


1 


15 


10.3999383 


4.586 


4.707 


1.7632 10^^ 


2.117 


8.006 


0.27733 


9.5000 10~ 


1 


20 


6.3947664 


4.879 


4.758 


1.9554 10^1 


2.001 


8,018 


0,34013 


9.5000 10~ 


1 


40 


2.4225710 


5.507 


4.856 


2,5643 10^1 


1.767 


8,043 


0,46558 


9,5000 10" 


1 


100 


1.1063569 


6.203 


4.935 


3.9656 10^1 


1.529 


8.075 


0,57258 


9.5179 10" 


1 


200 


0.9510749 


6.633 


4.950 


6.0802 10^^ 


1.356 


8.089 


0.53188 


9.4999 10~ 


1 


300 


0.9983394 


6.866 


4.946 


8.1006 10-^^ 


1.263 


8.096 


0.48495 


9.4997 10~ 


, 1 


400 


1.1107156 


7.027 


4.933 


1.0328 10^^ 


1.200 


8.101 


0.42339 


9.5000 10~ 


1 


600 


0.9687933 


7.246 


4.963 


1.1603 10^^ 


1.121 


8.109 


0.45154 


9.4998 10~ 


1 


Central Hydrogcn= or 


Scattering St 


allation mode. 


1 












5 " 


5.0000000 


2.830 


4.322 


1.3769 10^^ 


1.384 


7.452 


0.75400 


4.8940 10" 


■2 


7 


60.2205803 


3.758 


4.563 


1.3173 10^^ 


2.903 


8.048 


0.00000 


0.0000 




9 


32.9417891 


4.076 


4.603 


1.5812 10^^ 


2.758 


8.079 


0.00000 


0.0000 




12 


17.9675788 


4.441 


4.642 


2.0086 10^^ 


2.615 


8.109 


0.00000 


0.0000 




15 


12.6592264 


4.728 


4.665 


2.5143 10^^ 


2.517 


8.131 


0.00000 


0.0000 




20 


8.6058962 


5.038 


4.695 


3.1378 10^^ 


2.414 


8.151 


0.00000 


0.0000 




40 


4.3902484 


5.670 


4.740 


5,2772 10^1 


2.218 


8,190 


0,00000 


0,0000 




100 


2.9876710 


6.358 


4.609 


2,1264 10^^ 


1.978 


8,225 


0,00000 


0.0000 




200 


2.3571500 


6.730 


4.639 


2.8487 10^^ 


1.834 


8.240 


0.00000 


0.0000 




300 


2.1357779 


6.942 


4.651 


3.4292 10^^ 


1.747 


8.248 


0.00000 


0.0000 




400 


2.0465587 


7.106 


4.637 


4.4207 10^^ 


1.691 


8.255 


0.00000 


0.0000 




600 


2.0407706 


7.362 


4.810 


2.6816 10^^ 


1.595 


8.265 


0.00000 


0.0000 




^ Track Stalled 


















Table A4. AC and SC mechanisms active; 






1 n— 38 

lu 


cm 






M/Mq 


Age(Myr) 


Log(L/Lg) 


Log(T^j.j.) 


1 R.* (cm) 




Log(Tc 


) Xc 






Hayashi M 


inimum Lumit 


rosity Model 
















5 


0.0732837 


2.865 


3.859 


1.2047 10^^ 


-0.975 


6.618 


- 


- 




7 


0.0207489 


3.115 


3.730 


2.9142 10-*^^ 


-2.650 


6.130 








9 


0.0177888 


3.449 


3.724 


4.3936 10^^ 


-3.036 


6.069 








12 


0.0153947 


3.814 


3.717 


6.9172 10^^ 


-3,536 


5,982 








15 


0.0144330 


4.089 


3.713 


9.6710 10^2 


-3,842 


5,936 








20 


0.0124400 


4.410 


3.707 


1.4378 10^^ 


-4,283 


5.859 








40 


0.0082365 


5.094 


3.699 


3.2837 10^^ 


-5.150 


5.707 








100 


0.0043604 


5.805 


3.704 


7.2797 10^^ 


-5.820 


5.613 








200 


0.0030809 


6.284 


3.723 


1.1551 10^^ 


-6.077 


5.605 








400 


0.0016042 


6.745 


3.748 


1.7519 10^"* 


-6.234 


5.618 








600 


0.0004468 


6.998 


3.756 


2.2582 10^"* 


-6.386 


5.603 








Ln„o./L. 


= 95% 


















40 


20.7324568 


5.628 


4.754 


4.7132 10^^ 


1.800 


8.058 


0.10605 


9.5000 10" 


1 


100 


4.5980377 


6.279 


4.792 


8.3431 10^^ 


1.540 


8.080 


0.12656 


9.5133 10" 


1 


200 


3.5912405 


6.718 


4.748 


1.7011 10^^ 


1,375 


8,094 


0,14434 


9.5000 10" 


1 


400 


3.4325657 


7.085 


4.740 


2.6877 10^2 


1,241 


8,110 


0,11646 


9,4998 10" 


1 


600 


3.0620294 


7,308 


.1,6. '2 


.5.2214 10^^ 


1,18.' 


8,122 


0,0840.' 


9, ,^,000 10" 


1 




5 « 


5.0000000 


2.698 


3.732 


1.7904 10^2 


-2.567 


6.154 


0.75500 


0.0000 




7 " 


5.0000000 


3.191 


3.956 


1.1241 10^^ 


-1.067 


6.721 


0.75500 


1.5187 10" 


8 


9 " 


5.0000000 


3.578 


4.119 


8.2850 10^^ 


-0,652 


6,925 


0,75500 


3.4554 10" 


■7 


12 * 


5.0000000 


3.984 


4.309 


5.5088 10^^ 


-0,068 


7,190 


0,75500 


1,6556 10" 




15 " 


5.0000000 


4.274 


4.438 


4.2336 10^^ 


0.342 


7.376 


0.75500 


2.6014 10" 


4 


20 " 


5.0000000 


4.630 


4.581 


3.3130 10^1 


0.774 


7.578 


0.75440 


2.2035 10" 


3 


40 


21.0633259 


5.674 


4.735 


5.4184 10^^ 


2.218 


8.190 


0.00000 


0.0000 




100 


4.8774489 


6.311 


4.749 


1.0565 10^^ 


1.984 


8.220 


0.00000 


0.0000 




200 


3.8658341 


6.744 


4.619 


3.1638 10^^ 


1.833 


8.241 


0.00000 


0.0000 




400 


3.6641729 


7.105 


4.454 


1.0244 10^^ 


1.702 


8.256 


0.00000 


0.0000 




600 


4.1027411 


7.322 


4.418 


1.5577 10^3 


1.626 


8.265 


0.00000 


0.0000 





Track Stalled 



Table A5. AC mechanism only, for a neutralino mass mj(.=200GeV. 



M/M^ 


Age(Myr) 


Log(L,/LQ ) 




R*(cm) 


Log(pc) 


Log(Te) 










Starting P. 


aint 






















5 


0.0000000 


4.882 


3.642 


3.3438 


10l3 


-6.576 


4.810 


- 




- 




7 


0.0000000 


5.009 


3.647 


3.7899 


lo" 


-6.601 


4.901 


- 




- 




9 


0.0000000 


5.104 


3.649 


4.1733 


ioi3 


-6.615 


4.960 


- 




- 




12 


0.0000000 


5.212 


3.653 


4.6489 


ioi3 


-6.615 


5.025 


- 




- 




15 


0.0000000 


5.236 


3.659 


4.6423 


10l3 


-6.465 


5.119 


- 




- 




20 


0.0000000 


5.417 


3.660 


5.6860 


j^qIS 


-6.603 


5.129 


- 




- 




40 


0.0000000 


5.699 


3.673 


7.4259 


10l3 


-6.532 


5.264 


- 




- 




100 


0.0000000 


6.099 


3.698 


1.0491 


lO^'* 




5.417 










200 


0.0000000 


6.400 


3.716 


1.3640 


loi* 


-6.370 


5.508 


- 




- 




41)1) 
(lUl) 




- u-u 


:;.7r,2 


., ij-'^-- 


11)1-1 
11)11 


-n. 1^7 


.^..033 












= 50% 






















5 


0.0080226 


3.090 


3.716 


3.0198 


ioi= 


-3.244 


5.928 


- 




- 




7 


0.0123634 


3.216 


3.722 


3.3996 


10l2 


-3.128 


6.007 


- 




- 




9 


0.0139301 


3.452 


3.723 


4.4418 


10l2 


-3.139 


6.040 


- 




- 




12 


0.0140096 


3.872 


3.724 


7.1618 


10l2 


-3.235 


6.074 


- 




- 




15 


0.0116689 


4.174 


3.724 


1.0123 


10l3 


-3.451 


6.060 


- 




- 




20 


0.0108273 


4.537 


3.738 




jgl3 


-3.669 


6.059 


- 




- 




40 


0.0070071 


5.266 


3.778 


2.7769 


10l3 


-4.223 


6.017 


- 




- 




100 


0.0039907 


6.043 


3.975 


2.7385 


10^3 


-3.902 


6.255 










200 


0.0026829 


6.516 


4.059 


3.2085 


10l3 


-3.826 


6.355 


- 




- 




400 


0.0016259 


6.935 


4.111 


4.0936 


10l3 


-3.796 


6.431 










600 


0.0013353 


7.164 


4.129 


4.9076 


10l3 


-3.795 


6.467 












10% 






















5 


0.0509513 


2.806 


3.765 


1.7354 


10l2 


-1.431 


6.428 


- 




- 




7 


0.0406938 


3.361 


3.907 


1.7131 


10l2 


-1.271 


6.585 


- 




- 




9 


0.0315087 


3.723 


3.980 


1.8568 


10l2 


-1.318 


6.642 


- 




- 




12 


0.0257831 


4.113 


4.081 


1.8228 


10l2 


-1.248 


6.745 


- 




- 




15 


0.0207747 


4.390 


4.144 


1.8789 


10l2 


-1.240 


6.806 


- 




- 




20 


0.0182315 


4.720 


4.235 




10l2 


-1.131 


6.912 


- 




- 




40 


0.0115195 


5.392 


4.374 


2.0643 


10l2 


-1.133 


7.053 


- 




- 




100 


0.0063115 






3.1661 


loi^ 














200 


0.0044895 


6.553 


4.521 


3.9951 


10l2 


-1.346 


7.182 


- 




- 




400 


0.0032827 


6.958 


4.559 


5.3438 


10l2 


-1.387 


7.234 










600 


0.0029023 


7. 180 


4.574 


6.4304 


10l2 


-1.418 


7.259 










Hayashi M 


inimum Lumii 


losity Model 




















5 


0.0300699 


2.650 


3.739 


1.6362 


ioi2 


-2.059 


6.233 


- 




- 




7 


0.0182169 


3.115 


3.730 


2.9145 


ioi2 


-2.650 


6.130 


- 




- 




9 


0.0148365 


3.449 


3.724 


4.3933 


ioi2 


-3.037 


6.068 


- 




- 




12 


0.0122626 


3.814 


3.717 


6.9151 


ioi2 


-3.539 


5.981 


- 




- 




15 


0.0101071 


4.088 


3.713 


9.6724 


loi^ 


-3.855 


5.932 


- 




- 




20 


0.0091717 


4.409 


3.707 


1.4361 


ioi3 


-4.275 


5.861 


- 




- 




40 


0.0061004 


5.091 


3.700 


3.2531 


lo" 


-5.074 


5.731 


- 




- 




100 


0.0026849 


5.807 


3.703 


7.3077 


jq13 


-5.838 


5.608 










200 


0.0009983 


6.297 


3.721 


1.1832 


ioi4 


-6.124 


5.589 


- 




- 




400 


0.0000092 


6.761 


3.752 


1.7546 


ioi4 


-6.186 


5.634 










600 


0.0000092 


7.029 


3.786 


2.0349 


ioi4 


-6.002 


5.730 










Ln„c./L. 


= 95% 






















5 


0.6992481 


2.927 


4.452 


8.4513 


lO^" 


2.080 


7.673 


0.75283 


9 


.5000 10~ 




7 


0.4522332 


3.380 


4.551 


8.9905 


lolo 


2.108 


7.775 


0.75226 


9 


.5000 10~ 


1 


9 


0.3378005 


3.705 


4.623 


9.3935 


ioi° 


2.137 


7.850 


0.75159 


9 


.5000 10~ 


1 


12 


0.2151746 


4.071 


4.703 


9.9326 


loi" 


2.170 


7.931 


0.75200 


9 


.5000 10~ 


1 


15 


0.1683067 


4.342 


4.762 


1.0317 


loii 


2.199 


7.992 


0.75207 


9 


.5000 10^ 


1 


20 


0.1110554 


4.699 


4.829 


1.1438 


lo" 


2.201 


8.051 


0.75310 


9 


.5000 10~ 


1 


40 


0.0483624 


5.455 


4.921 


1.7885 


loii 


1.999 


8.100 


0.75440 


9 


.5000 10~ 


1 


100 


0.0261951 


6.161 


4.988 


2.9566 


lo" 


1.749 


8.139 


0.75480 


9 


.5002 10~ 


1 


200 


0.0193638 


6.599 


5.018 


4.2712 


loii 


1.586 


8.159 


0.75480 


9 


.5003 10~ 


1 


400 


0.0156433 


6.994 


5.033 


6.2776 


loii 


1.443 


8.177 


0.75490 


9 


.6414 10~ 


1 


600 


0.0142502 


7.208 


5.039 


7.8114 


loii 


1.361 


8.186 


0.75490 


9 


.6439 10~ 


1 


Central Hydrogen= or 


Scattering Stallation model 


















5 


56.1384003 


3.340 


4.505 


1.0644 


lo" 


3.200 


7.983 


0.00000 




0.0000 




7 


29.7313755 


3.734 


4.566 


1.2641 


loii 


2.924 


8.040 


0.00000 




0.0000 




9 


19.5580299 


4.053 


4.604 


1.5289 


lo" 


2.771 


8.076 


0.00000 




0.0000 




12 


13.7294256 


4.447 


4.640 


2.0425 


loii 


2.613 


8.110 


0.00000 




0.0000 




15 


10.7102756 


4.726 


4.665 


2.5130 


lo" 


2.518 


8.131 


0.00000 




0.0000 




20 


7.6979362 


5.040 


4.694 


3.1552 


loii 


2.415 


8.152 


0.00000 




0.0000 




40 


4.1505195 


5.668 


4.741 


5.2413 


loii 


2.222 


8.190 


0.00000 




0.0000 




100 


2.7010537 


6.327 


4.735 


1.1476 


10l2 


1.982 


8.221 


0.00000 




0.0000 




200 


2.2941418 


6.744 


4.547 


4.4148 


ioi2 


1.843 


8.242 


0.00000 




0.0000 




400 " 


1.9648883 


7.194 


4.342 


1.9067 


10l3 


1.191 


8.105 


0.16920 


9 


.7490 10~ 


1 


600 


2.0679154 


7.538 


5.111 


8.2005 


lo" 


1.583 


8.275 


0.00000 




0.0000 





°- Problems to reach end of H-burning 



